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Abstract. The energy, area and excess energy of a decorated vertex in a 2D foam are calculated. The 
general shape of the vertex and its decoration are described analytically by a reference pattern mapped 
by a parametric Moebius transformation. A single parameter of control allows to describe, in a common 
framework, different types of decorations, by liquid triangles or 3-sided bubbles, and other non-conventional 
cells. A solution is proposed to explain the stability threshold in the flower problem. 

PACS. 83.80.Iz Emulsions and foams - 82.70.Rr Aerosols and foams - 82.70.Kj Emulsions and suspensions 
- 68.03.Hj Gas-liquid and vacuum-liquid interfaces: Structure, measurements and simulations 



1 Introduction 

Since J. Plateau p], 2D foams have been extensively stud- 
ied, experimentally and theoretically, because they are 
simpler than three-dimensional systems [2J. 

In the dry model of 2D soap foams, the gas is assumed 
incompressible and the liquid fraction is assumed vanish- 
ing. Although the dry model correctly describes some as- 
pects of foam physics like energy minima, equilibrium con- 
figurations, etc 2 , in many other cases, the presence of the 
liquid needs to be taken into account, both theoretically 
[3] and to match experimental observations |4l5l6j . 

The more realistic model of foams, including a small 
liquid fraction at the vertices, is related to the dry model 
by Weaire's decoration theorem. The decoration theorem 
[7] state that, in 2D foams, the films connecting to a 
Plateau border with vanishing disjoining pressure, if con- 
tinued inside it, all intersect at a single point and satisfy 
the equilibrium conditions (fig. QJa). Conversely, it is al- 
ways possible to decorate a 3-fold vertex at equilibrium 
by a Plateau border, that is, to replace the vertex by a 
(small) triangle of liquid without changing the geometry 
outside the triangle and still satisfying equilibrium. 

The star-triangle equivalence [3] extends the decora- 
tion theorem to general cellular systems where the surface 
tensions have arbitrary values; this includes 3-sided bub- 
ble in standard foams (constant uniform surface tension) 
or in non standard foams (where the surface tension takes 
different values on different films [8 9J), or Plateau bor- 
ders with non vanishing disjoining pressure. In this con- 
text, bubbles, drops and Plateau borders are treated on a 
common footing as cells. 
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The decoration theorem validates the ideal dry model 
for slightly wet foams as far as equilibrium is concerned [2] ; 
the star-triangle equivalence permits, in principle, to take 
away 3-sided bubbles in searching equilibrium configura- 
tions. However, the foam energy, the area of the bubbles, 
coarsening and most mechanical properties of the foam 
change if the decoration is switched on or off. 

In dry foams, since the energy is the product of surface 
tension by film length, minimisation of line-length at fixed 
bubbles area completely determines the 2D foams equilib- 
rium structure. The energy of a progressively strained dry 
foam increases until two threefold vertices approach one 
another and undergo a Tl neighbour exchange [10111112] . 
The Tl quickly reduces the energy by a finite amount: 
the energy difference between the configuration with an 
unstable fourfold vertex and the one with a pair of three- 
fold vertices linked by a film. That is why the location and 
statistics of Tl events determine the inelastic properties 
of the foam and how it releases energy under strain or 
quasi-static flow. 

If the foam is not dry, the presence of Plateau bor- 
ders reduces the film length, and thus favours the switch 
compared to a dry foam [13]. At the switch, two triangu- 
lar Plateau borders meet, merge [14115] and then split to 
a new configuration involving again two distinct Plateau 
borders. So, in terms of energy furnished to the foam, the 
presence of Plateau borders reduces the barrier one needs 
to overcome to trigger a Tl. Furthermore, it is well known 
that the yield stress of a foam decreases with the liquid 
fraction [I6ll7ll8ll9j . 

By similar considerations, when 3-sided bubbles are 
inserted at some vertices of the foam, one can expect an 
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increase of the Tl energy barrier and, therefore, a reduc- 
tion of the shear plasticity due to edge flips. 

In this paper, we calculate the energy of a wet-deco- 
rated foam starting from the energy of the dry model. 

Recently, Teixeira and Fortes [20] gave the equations 
describing the exact geometry of a general Plateau bor- 
der with zero disjoining pressure: the surface tension of 
the films and of the liquid-gas interfaces are related by 
Tfiim = 27bordcr; they calculated the excess energy among 
other quantities. The excess energy is the energy difference 
between a decorated and a dry vertex. 

Using the invariance of 2D foams under Moebius trans- 
formations [21], we generalise this problem to the cases 
where the surface tension of the films and that of the dec- 
orating triangles are arbitrary, giving the Plateau border 
and 3-sided bubble problems a unified description. We cal- 
culate the excess energy E and all the geometrical quanti- 
ties as functions of four parameters which characterise the 
cluster size and shape. Normalising the excess energy by 
the square root of the triangle area gives a scale invariant 
quantity ep: the relative excess energy. 

In realistic foams, the relative excess energy of 3-cells 
is found to be approximately linear as a function of the 
triangle area. The slope of ep depends on the form of the 
decorated films. According to Lewis' law [22], the average 
cell area increases linearly with the number of sides or 
neighbours: (A n ) cx n. So the triangle area A3 is small on 
average. As we'll see, the slope is proportional to the sum 
of the squared curvatures of the films meeting at the deco- 
rated vertex. Being zero in the 3-fold symmetric case, the 
slope of ep measures the deviation from perfect symmetry. 

In the final part, we apply our analytical results to 
the flower problem. In the experiments, a bubble spon- 
taneously gets out [5] before the critical area is reached 
where the dry model predicts a spontaneous symmetry 
reduction [23]. Following the experimentalists' suggestion 
that the ejection might be due to the presence of liquid, 
we solve the problem by including liquid triangles around 
the vertices and calculating a threshold area for the cen- 
tral bubble at which two vertices get into contact. Our 
conjecture is that this contact triggers the ejection. 

The paper is organised as follows: Section 2 recalls the 
equilibrium equations for 2D foams with non-constant sur- 
face tensions and the star-triangle equivalence. Section 3.1 
defines the reference structure which is then used, in sec. 
3.2 and 3.3, to describe general 3-cells, including Plateau 
borders, 3-sided bubbles and similar patterns. Section 4 
contains the calculation of the main quantities: the en- 
ergy gap between bare and decorated structures, the area 
of the 3-cell and the relative excess energy. In section 4.4, 
the relative excess energy is studied as a function of the 
pressure of the four bubbles involved. Section 5 contains 
the series expansion of the relative excess energy with re- 
spect to the 3-cell area when the star films are fixed. Fi- 
nally, in Section 6, we apply our analytical results to the 
flower problem. 



2 Equilibrium of 2D foams 

To describe various types of 2D foams, we consider here a 
general 2D cellular system with surface tensions verifying 
the following equilibrium properties. 

The cells are separated by interfaces obeying Laplace- 
Young's law: 

AP + 7 fc = 0. (1) 

AP = P<l — P\ is the pressure difference across the in- 
terface, 7 the surface tension and k the curvature of the 
interface |2J . In static conditions and in absence of applied 
field, the pressure inside each cell is constant so that the 
interfaces are arcs of circlefQ. 

At equilibrium, the interfaces meet three by three at 
vertices in a manner satisfying Plateau's laws [112] : 

e? ,-,», (2) 

E- =1 7A=0, (3) 

where tj is the unit vector tangent to the interface at the 
vertex and where the surface tension, jj (jj > 0), differs 
from interface to interface }3I8I9|24| . 

Equations (J2J and ([3|) imply that the symmetry group 
of 2D cellular systems at equilibrium is the group of homo- 
graphies, or Moebius maps, generated by Euclidean sim- 
ilarities and inversion |24|25j . Some definitions and basic 
properties of Moebius maps are recalled in appendix [5] 

Equations {2j and ([3|) also imply that the centres of 
curvature of three circular edges meeting at a vertex are 
aligned [24] . 

This common framework describes several situations. 
A cell containing liquid is a Plateau border. A cell filled 
with gas is commonly called a bubble. The sides of a liquid 
cell are liquid-gas interfaces, with liquid-gas surface ten- 
sion 7^ g . The edges separating (gaseous) bubbles are liquid 
films, containing a negligible amount of liquid compared 
to Plateau borders. Films have a specific surface tension 
781m- In standard dry equilibrated 2D foams, the surface 
tension is constant over the entire foam (7^ = 7 = 7fiim)- 
In non-standard dry 2D foam, the surface tension may 
vary from film to film (7^ jji ) |8l9j . Slightly wet equili- 
brated 2D foams contain bubbles and a small amount of 
liquid mainly confined in Plateau borders forming concave 
triangles around the bubble vertices. 

The energy of a 2D cellular cluster is the sum, over the 
set of edges j, of the edge length, lj, weighted by surface 
(line in 2D) tension: 

£ = Ev,- (4) 

3 

2.1 Star-Triangle Equivalence 

A triangle is a three-sided equilibrated cell. Each vertex 
is the end of exactly 3 curved edges; two of them are part 

1 Even if we admit surface tension varying from edge to edge, 
we always assume that it remains constant inside every edge. 
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of the triangle boundary but the third one is external, 
outside the triangle. If the foam is viewed as a graph, this 
edge connects the triangle to the rest of the foam. Each 
triangle has three such connecting edges, or bonds, one at 
each vertex. 

A star is an equilibrated figure formed by a vertex and 
the three edges meeting at it. 

At equilibrium the pressure, or the area, and the sur- 
face tensions fix the triangle geometry. The star-triangle 
equivalence states the following [3J: 

— For any triangle, the three external bonds, if extended 
inside the triangle with the same curvature, all inter- 
sect at a common point, that forms a star with the 
extended edges. 

— Conversely, given any star, there is a continuous range 
of areas extending down to 0, and an open set of ten- 
sion values, such that, for any prescribed values in 
these ranges, there is a triangle with one vertex on 
each branch of the star that forms an equilibrated fig- 
ure once the portion of the star inside the triangle has 
been removed. Briefly said, the triangle decorates the 
star. 





b) 

Fig. 1. Star-triangle equivalence in the case of a Plateau bor- 
der (a) and a three-sided bubble (b). The triangle (solid lines) 
can be replaced by the star (dotted) and, conversely, the star 
can be decorated by the triangle, in a way preserving equi- 
librium. The connecting films (in grey) around the triangle 
continue, as dotted arcs of circles, inside the triangle to meet 
at the equilibrated point vq. 



According to this equivalence, the star and triangle can re- 
place one another without any modification in the rest of 
the foam, and in a way preserving all the equilibrium con- 
ditions. To be precise, both figures have the same contact 
points with the outside foam, namely the triangle vertices; 
in the substitution, the star branches are limited to the 
segments inside the triangle. The external parts are kept 
unchanged. The star-triangle equivalence was proved in 
full generality for cellular systems with multiple surface 
tensions [3j. 

When we consider a triangle and its associated star as 
a whole, a star+triangle, we will call it 3-cell. We need the 
complete figure to evaluate the excess energy, for example. 



But it should stay clear that, physically, the star and tri- 
angle cannot be both present at the same time. The only 
exceptions will be found in sec. 13.31 

The star-triangle equivalence contains Weaire's deco- 
ration theorem [7] by triangular Plateau borders as a par- 
ticular case (fig.[I]-a): when the surface tension of the trian- 
gle sides is that of a liquid-gas interface, j e = 7 tr i = 7^, 
while the star has the film surface tension, 7, = 7 sta = 
7fii m ~ 27^ 9 , at zero disjoining pressure. 

For a standard dry foam, star-triangle equivalence in- 
volves a star and a 3-sided bubble with the same tension 
everywhere: 7, = ~/ e = jshn. 

In 3 dimensions, a decoration-bubble theorem holds 
only for spherical foams [26] . When the films are not nec- 
essarily spherical, Teixeira and Fortes proposed a modi- 
fied, approximate, version involving both line and surface 
tensions l27l. 



3 Conformal description of 3-cells 

In this section the description of the 3-cell (star+triangle) 
is given. First, we construct a reference 3-cell depending 
on two parameters. Then, a general 3-cell is obtained by 
a Moebius transformation. 

To fix notations (see fig[T|), let ip tT1 , of radius r trl , be 
the edges of the triangle (decorating edges), and ip sta , of 
radius r sta , the edges of the star (decorated or internal 
edges). When the star, with vertex v , is decorated by 
a triangle, vq and ip sta are virtual. Conversely, when the 
star is not decorated, the </? trl are virtual. We will use the 
convention of signed angles and radii, in a way such that 
the arc length is always positive [7]. 

In this paper, we allow only two different values for 
7: the surface tension on the triangle boundary sides is 
7 e = 7t r i, and that of the external and internal films, on 
the star edges, is 7, = 7 s t a - These tensions are related by 
equation ([2]), applied to a vertex of the triangle: 



7i = 2 7 e cosa. 



(5) 



The contact angle a is the angle between (p sta and ip trl . 

If the thickness of the films, h, is not negligible, the 
angle a is related to the disjoining pressure II [819126) : 



n 



li, 

— (cos a 



1) 



The excess energy is the energy gained by decorating 
the vertex. It is the difference of the internal and triangle 
edge lengths, weighted by the surface tensions [20] : 



E — -^triangle 
L tri 

li 



£star = Tei' 1 ' 1 ~ li^ 



2 cos a 



(6) 



The relative excess energy is defined as the ratio of the 
excess energy E over the square root of the triangle area: 

(7) 



~ A 1/2 

triangle 

From now on, we will set 7, = 1; the surface tension 
of the internal films is our tension unit. 



4 



M. Mancini, C. Oguey: Decorated vertices with 3-edged cells in 2D foams: exact solutions and properties. 



3.1 The reference 3-cell 

The reference 3-cell C(q,f3), of parameters (3 and q, is the 
completely symmetric star+triangle defined in figure O 




Fig. 2. The reference, symmetric, 3-cell with parameters j3 
and q. The star (internal films) is represented by dashed lines. 
VC is a radius of the left side. Solving the triangle VOC shows 
that the angle at C is indeed (5 = a — 7r/6. The case of a bubble, 
ft > 0, is on top; the Plateau border case, (3 < 0, is below. 



The parameter q > is the chord length of the curved 
sides of the reference triangle. It specifies the triangle size. 
The parameter (3 is the angle defined by 

(3 = a — 7r/6 . 

As conformal means angle preserving, the contact angle 
a, and so (3, are preserved by the Moebius transformation 
applied to, or from, the reference cell. Choosing (3, rather 
than a, helps in treating Plateau borders and 3-sided bub- 
bles in a common way. In figure (JTJ-a), (3 is equal to 7r/6; 
in fig. |T]-b), (3 = — 7r/6. In the reference cell, the angle 
subtended by the triangle sides equals |2/3|; this property 
is not generally preserved by conformal transformations. 

When (3 = 7r/3, ji/je — and the sides of the triangle 
C(q,f3) form a circle. This case describes a 3 cells parti- 
tion of a disc surrounded by a rigid membrane. So, the 
interpretation of the triangle as a bubble or a liquid drop 
is limited to the values (3 £ [— 7r/6, 7t/3]. 

The edge curves are parametrised by 



ipf a (h) = e lS > h , with ti e [o,q/Va], 



qe 



2sin(/3) VV3 



— p= cos(/3 + 7r/6) — e 



(8) 
, (9) 



with to E [—1,1]: j = 1,2,3 indexes the edges and 9j = 
2nj/3. 

The curve © depends continuously on (3 € [— 7t/6, 7t/3] . 
Its curvature vanishes, and changes sign, at $ = 0. The 
equilateral triangle is curved, convex, at positive /?; curved, 
concave, at negative /3; and it has straight edges at (3 = 0. 
At f3 = the parametrisation is given by the limit of ([9]) 
as (3 -> 0: 

^f{t 2 ) = -qe^ l -(i- 1 ' 2 +it 2 ) . (10) 

With © , the excess energy of the reference cell is 

3/3 



2 sin (3 cos (/? + f) 



V3 



q- 



(ii) 



The triangle area is the sum of the area, , of the 
rectilinear equilateral triangle based on the same vertices, 
and Au, the signed area of the three lenses around the 
straight triangle: 



. , , , . / V3 3 2(3 - sin(2/3) \ , 
Mq)=A A+Au = (- + - \f. 



(12) 



Finally, the reference relative excess energy is obtained 
by dividing the excess energy by the area square root: 



EM 
Ap{q) 1/2 



'3/3 -2^3 cos(/3+£) sin/3 

^ sis °<« J ^WTT) • < 13 > 

Being dimensionless, does not depend on q. In the par- 
ticular cases of a Plateau border and of a 3-sided bubble, 
the respective are 



, ; - - v V3-tt/2~ -0.401565, 

4 /a = ^/2(7r-V3) - 1.67901. 

As a final remark on the reference, notice that the star 
is fixed, made of straight edges meeting at 27r/3, the stan- 
dard Plateau angles. In particular, the star is independent 
of the parameter (3 and it is affected by q only regarding 
the length at which its edges are cut. This indicates why, 
in all subsequent formulae, the quantities concerning only 
the star do not depend on /3, nor on q except for side 
length. 



3.2 The general 3-cell 

Applying a suitable Moebius transformation to the refer- 
ence 3-cell produces a general 3-cell at equilibrium (the 
star and triangle both satisfy eqs. (P), (f5]) and ([3])). 

Given a star, the decorating triangle is entirely speci- 
fied by its area A and contact angle a; those will be related 
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to the parameters q and (3 of the reference 3-cell. Now, our 
transformation must map the (fixed) reference star onto 
a general star of any possible shape (at equilibrium, of 
course). But, in the star, centred at vertex v°, the three 
(internal) films i — 1,2,3 have curvature fc^. Because of 
equation ([3]), only two of the curvatures are independent. 
So, removing the Euclidean degrees of freedom, we need 
only a one-complex parameter Moebius transformation to 
reach all possible star+triangle figures. 

A general Moebius transformation /, given by (|44|). has 
six real parameters, much beyond our needs. We require 
that the transformation i) preserves orientation and ii) 
maps the interior of the reference 3-cell onto the interior of 
the general 3-cell. In addition, we can fix a length scale, the 
star vertex position, and one of the film tangents there (by 
conformality, this fixes all the tangents at v°). Choosing 
the origin at v°, this means /(0) = and f'(0) > 0. Then 
the Moebius transformation can be written in the form 

/ s (z)-(l-| S | 2 )-^, with | S | 2 = S 5<1, (14) 
1 — s z 

depending on the single complex parameter s. The result 
is a 4 (real) parameters star+triangle figure: C(s, q, (3) = 
f s (C(q, (3)). We designate the transformed quantities by a 
tilde (appendix |A|) . An example is given in figure [3l 




Fig. 3. a) The reference 3-cell C(q,/3), (3 = tt/6 (in black) 
and /3 = —n/6 (in grey), b) General 3-cell C(s,q,(3), image 
of a) by f s . The figures are calculated for q = y/2 and s — 
0.2 exp (— i-|). The circumscribed circle in b) is given by (|15[1 . 



The points z — and s are fixed points of f s . The com- 
plex function f s maps the reference circumscribed circle, 
of centre Co = «o = and radius r$ = qj V3 (fig. [3]), onto 
the circle of parameters (eq. (j45|l): 

( r ~, ( 2 f(i -N 2 ) V3(i-H 2 ) \ 

(Co ' ro) = \ q 3^W 9 |3-^| S p| ) ■ (15) 

The 3-cell is contained in the disc (Co,r ) for \s\q < 
V3, while it is outside0 the disc for \s \ q > 

As the equilibrium equations and are left in- 
variant by f s , the image of the reference, either as a star 
or as a triangle, is in mechanical equilibrium. 

In order to get physically meaningful patterns, q must 
be bounded: q < q max (s,/3). An expression for q max (s,/3) 

2 This peculiar situation is possible only when the triangle 
is curved and has two concave sides. 



will be derived in section 14.11 Here, we just give its ori- 
gin and meaning. At q = g max , the reference boundary, 
dC(q,f3), meets 1/s, the pole of / s ; this condition de- 
fines q max (s,/3). When q < q max , the reference triangle 
is mapped to a bounded triangle but, when q > g max , the 
reference interior is mapped to the exterior of the triangle 
boundary, of infinite area, violating our requirement ii), 
above. 

The divergence at q — > g max is illustrated in figure 
IU First, recall that any equilibrated star vertex has a 




Fig. 4. Two different 3-cells C(s,q,/3) calculated for f3 = n/6, 
q and q , with q > q. The continuations of the internal films 
is plotted in grey (green on-line). 

conjugate v* where the circles supporting the star edges 
also meet. The position of v* is the mirror reflection of v 
through the line of the curvature centres [24l26j . For the 
reference C(q,j3), the conjugate of vq — is Vq — oo; in 
C(s,q,/3), Vq is given by 

v* =f s (o O ) = -s(\s\- 2 -l). (16) 

By increasing q alone, the triangle vertices move along 
the star edges, towards infinity in the reference, and so, 
towards Vq under the mapping f s (figure I3J . When q > 
<Zmax, the triangle area is infinite and all the edge radii are 
negative. 

In real foams, where the star edges end at neighbour 
vertices, this divergence is not reached. The route to it 
is deviated by topological changes. Indeed, suppose that 
the area of a 3-sided bubble is increased progressively, for 
example by injecting gas; much before one sees its area di- 
verging, topological changes will occur, either transform- 
ing the triangle into a higher polygon, or breaking films. 



3.3 Extension to special 3-cells 

Although the 3-cell is well defined only for values —tt/6 < 
[3 < 7r/3, the equations can be extended, for particular 
values of the parameters, to cover other physical or math- 
ematical situations. These cases are special because the 
star and triangle are taken here altogether to form a com- 
plete equilibrated cluster, without any outside foam. With 
contracting tensions, this is only possible for obtuse con- 
tact angles, that is, (3 > tt/3. 



6 



M. Mancini, C. Oguey: Decorated vertices with 3-edged cells in 2D foams: exact solutions and properties. 



3.3.1 Case (3 — ir/3: triple partition of the disc 

As noticed in sec. l3.lt the reference C(q, (3 — 7r/3) is a disc. 
f s maps it to a disc divided into three bubbles by films 
meeting the outer boundary orthogonalljQ. Cafiete and 
Ritore [28] proved that this type of graph is the unique 
least-perimeter way of partitioning the unit disk into three 
regions of prescribed areas. 




Fig. 5. For q = V3 and f3 = tt/3, varying a, the star describes 
all the least-perimeter partitions of the unit disk into three 
regions of prescribed areas. 



3.3.2 Case (3 = it/2: three bubble cluster 

When (3 = it/2, considering ip sta and ip tT1 as real films, the 
graphics of the reference 3-cell is equivalent to a symmetric 
three bubble cluster (figure Oa). The absolute value of 
IT3|) gives its relative energy [E/yfA): 

£ 3 = |e°/ 2 l = (4\/3 + 67r) 1 / 2 ~ 5.07718. 

Applying (fT4"l) to the reference produces a general three 
bubble cluster with constant surface tension (figure [6]-b) . 





b) 

Fig. 6. a) At (3 = 7r/2, the reference describes a cluster of 3 
identical bubbles, b) Three-bubble cluster transformed by f s . 



4 Area and energy of a general 3-cell 

To calculate the relative excess energy of a 3-cell C(s, q, (3), 
we need to calculate the star and triangle perimeters, its 
area and the edge curvatures. 



3 a = 7r/2 is the angle formed by the soap films with a 
smooth rigid wall or membrane. 



Let us start with the edge radii. Using equations (|45| 
and PrJ)) . they are, for j — 1, 2, 3, 



~sta _ 

3 2|s|sin( 



ft." 
j 



-9 S ) ' 

6 (1 - |s| 2 ) gsin/3 



4|3 2 sin/3 — qse 19 ? cos(/3 - 



f)l 2 -3<? 2 H 2 



(17) 
, (18) 



where S — arg(s). As anticipated (sec. 13. ip . the star radii 
f sta = f sta (s) don't depend on q. The radii inverse, that 
is, the curvatures, verify the equilibrium conditions ([3]), 
and the following equations: 



1 



3 

j 



= 6 



. (r ~sta )2 - V 1_ 1S | 

, 1 8|s| 2 cos 2 (f +13) 



(19) 



3 



1 vj ' 



4(|s|V cos(f -/3)-3sin/3) ; 
3(1- |s| 2 )2 g 2 



(20) 



Thus, the sum of the squared curvatures, in the star or 
in the triangle, doesn't depend on the argument of s. The 
only relevant parameter left for the star sum in eq. ()19|1 is 
\s\. 



4.1 Perimeters 

The value of the internal perimeter L sta (s, q) follows from 
a straightforward calculation detailed in appendix [Bj 



3 3 
L st \s,q) = J2 IT = E f f 



(21) 



where 



wf a (s,g) = 2arctan 



3-^ 2 q\s\ -cos(9j -6 S ) 



2 arctan 



cos(8j 



sin(9j — 9 a 



(22) 



is the (signed) angle subtended by the internal film <^ ta 

(fig. [7]). Once more, the internal perimeter L sta doesn't 
depend on [3. 

Analogously, the triangle perimeter is given by: 

3 3 

Lf(s, q) = J2 If = E r T( s > 9, 13) wf («, q, (3) (23) 

3=1 3=1 

where w] rl (s, q, (3) is the angle subtended by the triangle 
edge (pj 1 , defined by (|55]1 . The angles wj n verify the equa- 
tion 



3=1 



(24) 
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Fig. 7. Star and triangle for /3 = — tv/10. Summing the internal 
angles of the (solid) grey triangle gives (I24[) . 

consequence of the fact that the angles of the (rectilinear) 
triangle sum up to ir (see figure [7]). 

The limit value q max (s, (3) (sec. !3.2|) . where the bound- 
ary of the reference triangle meets the pole 1/s of f s , can 
be derived from (fTSl) . Indeed, fixing all the parameters (s 



(25) 



and (3) but q, the sign of kj = l/fj rl changes when the cir- 
cle supporting <^ rl meets the pole 1/s. For each j — 1, 2, 3, 
this occurs at the positive solution q — qj of 

k 1 = l/rf i (s,q,(3) = 0. 

The explicit expression of qj is: 



V3 [3 - 4 cos 2 ifi + f ) sm 2 (9j - 9 S )] 
2|s|cos(/3-f) 
\/3cos(/3+ f)cos(^ -9 S ) 
|s|cos(/3-£) ■ 



1/2 



(26) 



When g increases, the pole is first met by circle por- 
tions outside the triangle. The bound q = q max (s,/3) cor- 
responds to the encounter with the triangle boundary in 
the last of the three circles: 



i?max(s,/3) 



max qa 

3=1,2,3 



(27) 



For 



Q > Qn 



the interior of the 3-cell is no more 



is no longer verified, as illus- 



bounded and relation 
trated on fig. [4] 

4.2 Area 



The area of the 3-cell C(s,q,(3) is the sum of four parts: 
the rectilinear triangle based on the triangle vertices and 
the (signed) area of the three lenses around: 



Ap(s,q) = A A (s,q) + A u (s,q,/3) 



with 



A A = 2f t 1 ri f* ri sin(^) 



A 



sin(^) 



Mf -/3+f)| 



E=i(T) 2 (£i-sm^)- 



(28) 

(29) 
(30) 



In these equations, uij stands for <Dj rl . Although (3 appears 
in the RHS of (f2"9"| . it is clear that the area of the triangle 
A a doesn't depend on it. Indeed, the triangle A is built 
only on the vertices, independent of a, and so (3 (they are 
determined by q in the reference and then mapped by / s ). 



4.3 Relative excess energy 



Using (|2"1"T) . (|2"3"j) and ([25]). we obtain the relative excess 
energy for a general 3-cell C(s, q, (3): 



ep(s,q) 



Lf(s,q)/(2 cos(/? + tt/6)) 



,9) 



(31) 



Before we draw the plots, let us make a change of vari- 
ables to take advantage of the dilation invariance of the 
relative excess energy, already noticed. The new parametri- 
sation is defined as follows: 



p =<z(i-N 2 ) 

6 S = arg(s), 



with 
with 



P> 0, 

o < v < fj(P, 



(32) 



where fj(f3, 9 S ) — \s\q max (s, (3) is independent of \s\ by ([26]). 
The fourth parameter, (3, is unchanged. The parameter p 
corresponds to a magnification of the 3-cell, so the rela- 
tive excess energy doesn't depend on it. Moreover, from 
the symmetry of the reference 3-cell under the group D3 
generated by mirrors at 7r/3, the relative excess energy 
verifies: 

e/}(r),0 8 ) =6/3(77, 2tt/3 ± 8 ). (33) 




Fig. 8. Plot of the relative excess energy of a 3-sided bubble 
(/3 = 7r/6) as a function of 77 = q\s\ for some values of < 
@s < 7r/3. The thick solid line corresponds to 6 S = where the 
system has a mirror symmetry. The dot-dashed lines point out 
the limit values of calculated at fj. 

Figures [5] and O show the plots of 6/3(77, 9 S ) for a 3-sided 
bubble and a Plateau border, respectively. In both cases, 
the relative excess energy is minimal for 77 = 0, i.e. in the 
symmetrical reference configuration. 

In the 3-sided bubble case, for low values of 77, 6^/6 
increases with 9 S while, for greater 77, it decreases with 9 S ; 
the change, (d/d9 s )e(w/6,r],9 s ) = 0, occurs at 77 values 
close to 77(^/6, 0). At 77 = 77(71/6, 9 S ), when the area of the 
3-sided bubbles diverges (sec. 13. 2p . € 7r /e(fj(n/6,9 s ),9 s ) is 
constant as a function of 9 S , meaning that all the curves 
in fig. [8] end at the same ordinate. 

In the Plateau border case {(3 — —tt/6), e^/g is a non 
increasing function of 9 S : (d/d9 s )e^ 7l /6 < 0. For any 9 S , 
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<=_ T/6 =-VV3-ff/2— 0.4017 
£^/»(?,0)=0 




Fig. 9. The relative excess energy of a Plateau border (/3 = 
— 7r/6) as a function olr) = q\s\ for some values of < 9 S < 7r/3. 
The solid line corresponds to # s = 0, where the border has a 
mirror symmetry. 



at r] = r?(0), the triangle vertices stay on a straight line 
and the area of the rectilinear triangle is zero: Aa = 0. 
This implies that at i] = fj(0), the derivative (d/drjje-n/e 
has a singularity, corresponding to the cuspidal points 
of e- v /e(rj, 9 S ) in figure El Again, in this case, at rj = 
fj(— 7r/6, 9 s ) and S 7^ 0, when the area of the Plateau bor- 
der diverges (sec. I3.2[) . the relative excess energy doesn't 
depend on 9 S : e_ n /Q(fj, 6 S ) — 1.772 . . . (outside the figure). 



4.4 Pressure as coordinates 



Solving equations (|18[) and inserting into the formulae of 
sec. 14. Hl4~3l yields the relative excess energy as a function 
of the side radii f* rl , j = 1,2,3. By Laplace's equation 
([TJ , the curvatures are proportional to the pressure differ- 
ences. If Po denotes the pressure in the 3-cell and Pj the 
pressure in the neighbours j = 1,2,3, we can analyse ep 
as a function of the pressures. 

By the dilation invariance of e@, we can set the ab- 
solute value of one of the sides radii to 1 without loss of 
generality. So, fixing P — and r\" — sign(/3) , flj) implies 



Pi = 
P2 = 
P3 = 



sign(/3) 



2 cos(/3 + tt/6) ' 
1 

2f!p cos(/3 + tt/6) 
1 

2ff cos(/3 + tt/6) 



(34) 



In figures fTUl and [TT1 f/3 = 7r/6 and (3 = —tt/6, respec- 
tively) , we have drawn the parametric plots of the relative 
excess energy as a function of P3 for different values of P2 ■ 

In figure [TU1 the limit of e 7r / 6 (P 2j P3) when P 3 —> —00 
corresponds to configurations of a non-symmetrical cell 
when the area of the 3-sided bubble goes to zero. P2 = 
P3 — > — 00 gives the configuration "rat", whereas in all 
the other cases, P3 — > —00 at fixed P2, the bubble goes to 
the configuration "cat" (fig. [12]). 




"vr/6 




50 P3 100 

Fig. 10. Standard bubble case, /3 = 7r/6 and Pi = —1. Plot of 
e -n-/6 versus the pressure P3 for different values of P%. 




30 A 35 



Fig. 11. Plateau border, = —tt/6 and Pi = —1/2. e_^-/6 is 
plotted with respect to P3 for different values of P%. 



In the Plateau border case (figure ITT]) . when P2 = 
P3 — > +00, the border goes to a configuration "fox" (fig. 
T^j) . where one of the edges is straight. In the other cases, 
£ -tt/6 goes to zero. Obviously, in all cases, bubble or bor- 
der, equal pressure (Pi = P 2 = P3) gives ep(Pi,Pi) — e® 3 . 



5 Decorating a fixed star 

In this section we consider an ideal dry 2D foam at equilib- 
rium and we ask how much the energy of the foam varies if 
we replace one vertex by a Plateau border or a triangular 
bubble. 

This question is important in foam rheology [1011 3 29J . 
If the foam is not dry, the film length is reduced by the 
presence of the Plateau borders (see figure Q2]) , favour- 
ing the switch. Preliminary calculations on small samples 
show that the presence of liquid Plateau borders reduces 
the Tl energy barrier, whereas the presence of triangu- 
lar bubbles leads to locally more stable quadrangles, thus 
increasing an effective Tl barrier. Qualitatively, as com- 
pared to a bare foam, the net result is a decrease of the 
yield strain and stress when the foam is decorated by 
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cat 





_2_(V3 + ^l-2)=4.765 



e n/6~ I ==r 

V 3 JT/2+ V3 

r fox_ ^ft-V3)-" 



Fig. 12. Limit configurations found when one or two of the 
pressure differences with the neighbours go to infinity. "Rat" 
and "cat" describe 3-sided bubbles, while "fox" appears for 
Plateau borders. 



Plateau borders, but an increase when the decoration is 
by triangular bubbles. 

More extensive investigations are in progress to get the 
full rheology. Here, we show that the energy contribution 
due to decoration can be put into a simple tractable form. 




Fig. 13. The distance between two vertices in dry a foam, ld ry , 
is larger than in a wet foam, l we t- 



Varying the parameter q and fixing s, we can vary the 
decorating triangle area, Ap[s, q), keeping constant the 
star radii f^ ta . As said in the introduction, the size of the 
3-cells in real foams will, in general, be smaller than the 
typical size of other bubbles. In this approximation, we 
can write the relative excess energy as a function of the 
triangle area. To do so, we expand L tri , L sta and Ap(s,q) 
with respect to q, near q = 0: 



P A (s,q) = {l-\s\ 2 )(l}*q + 
L^(s,q) = (l-\s\ 2 )(L?»ql 
Ap(s,q) = (l-\s\ 2 f(A 2 q 2 



:f| s |y 



t sta 
^3 



0{q% (35) 
s\ 2 q a )+0(q% (36) 
0(q 5 )- (37) 



The coefficients of q n , which depend only on /3, are given 
in appendix [Cj Using the expansion to second order in 
(|37p to explicit q as a function of Ap and replacing this in 
the expression (|3"Tj) of ep(s, q), we find: 



e /3 



(s,Ap) = e%+rpK 2 ta {\s\)Ap 



o(if ), 



(38) 



where e° and K 2 ta are given by (fT3|) and (fl9|) , respectively. 
So the excess energy defined by Teixeira and Fortes [20] 
is, for small areas: 



E 



(s,V) 



el A 1 ' 2 + r K 2 ta (\s\) Af + 0(A 2 ). (39) 



The factor Fp in (f38|) . ([39]) depends only on (3 and its 
expression is: 



i 



6A 3 2 /2 



t tri 
^3 



tri Aj 
2A 2 



2cos(/3+ f ) 



7- sta 
^1 



A4 

2A 2 



/- sta 
^3 



(40) 



Given a fixed star in a 2D dry foam, the relative excess 
energy produced by decoration is, for small areas, linear 
in A. The angular coefficient involves the numerical factor 
rjg, depending on the type of the decorating 3-cell, and 
the star sum of the square curvatures K 2 ta . 

In particular, the relative excess energy depends on 
\s\ only through the product of the triangle area and the 
mean square curvature of the star: 



60 ( S ,V) = e /3 (A,} ■ K 2 ta (\ S \), 



(41) 



K^ ta does not de- 



In terms of the coordinates (|32j) . Ap 
pend on p. The argument 9 S appears only in higher order 

terms of the expansions of ep. In |38|) 

3 /2 

proportional to K s { a cos(3 s ). 

Figure [14] shows the plot of ep versus ApK. 
case of a three-sided bubble. 



the A 3 J 2 term is 



^2 
'■sta 



the 



1.73 










1.72 










1 .71 












0.2 


^^^0.4 


0.6 




1 .69 










1.68 


r e /r/6 









Fig. 14. The relative excess energy versus A/3ifs ta for different 
values of S at [3 = 7r/6. Dashed, the linear approximation of 



«/3 



For AaK Bta < 1, the error is smaller than 0.6%. 



The positive quantity K 2 ta , depending only on \s\, has 
the dimension of an inverse area. For a symmetrical star, 
K 2 ta (0) = 0, and it increases with \s\. Let us remark that, 
by ©, we can write: 



(42) 



K 2 ta = 2[(k 1 ) 2 -k 1 k 2 

where for simplicity we call kj = l/f^ ta . This quantity 
appears when 2D dry foams are studied using projective 
geometry and local coordinates associated to any vertex 
[261. 



+ (M 2 

l/ffa 
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6 The flower problem 

The flower problem, introduced by Weaire et al |23j . con- 
sists of a dry cluster of n equal-sized bubbles surrounding 
a central bubble of area A c (fig. [T5Ta'). By energy min- 
imisation using Surface Evolver, the authors computed 
the minimal energy configurations as a function of A c . 
They found that, for n > 6, the symmetrical configura- 
tion is minimal when A c is bigger than a critical area A* . 
When A c < A* c the symmetry of the cluster is reduced 
and the energy landscape contains multiple, but energet- 
ically equal, minima (as shown in [5]). The degeneracy 
follows from the spontaneously broken n-fold symmetry. 
The value of A* corresponds to a pressure of the central 
bubble equal to the pressure outside the cluster. 

In Cox et al's experimental set-up [S] , the bubbles were 
trapped between a glass plate and the top of the liquid. 
The area of the central bubble was slowly decreased but 
the predicted symmetry breaking was not observed. Be- 
fore the critical area, at some A* c * > A*, a Tl instabil- 
ity occurred and a bubble was ejected from the cluster. 
The authors explain the ejection by the presence a small 
amount of liquid, which provokes the Tl, and by some 
effects due to the third dimension. 

Considering only the effect of the liquid, we estimate 
its influence on the occurrence of topological changes and 
the petal ejection. In the dry 2D model, when a film length 
shrinks to zero, the unstable 4-fold vertex relaxes to a 
lower energy configuration where an edge joins two 3-fold 
vertices. 

Now, with a little liquid at the vertices, the film length 
is reduced by the presence of the Plateau borders, so that 
the topological change is likely to occur earlier. In the 
flower, because of the Plateau borders, the film length 
may vanish before the critical area A* c is reached. This is 
precisely confirmed by our calculations. 

We have decorated the flower vertices by Plateau bor- 
ders at constant pressure (fig. [TSl-b). The borders don't 
have the same area. Moreover by the presence of a liquid 
reservoir in the experimental set-up [5] , the liquid fraction 
is not fixed. Then, we calculated the value of the central 
area, A\i m , at which Plateau borders touch each other (fig. 
QjjJ-c), when the ejection is conjectured occur. 

As we didn't know the experimental liquid pressure, 
we chose the value, of the pressure in the Plateau bor- 
ders, that agreed best with the experimental points. This 
value is Po/li = — 0.390 ± 0.015 m _1 , provided the exter- 
nal pressure is set to zero. 




Fig. 15. Flower problem, a) The flower cluster without liquid, 
b) The vertices are decorated by Plateau borders, c) The limit 
where some Plateau borders touch each other. 




6 8 10 12 14 16 18 20 



Fig. 16. Critical area of the central bubble, normalised by 
the petal area, versus the number of petals n. Our calcula- 
tions are compared to the experimental data from [jjj. Squares: 
the calculated values at which two Plateau borders touch each 
other when the central area decreases at constant liquid pres- 
sure (.Aiim); these values are close to the solid curve of equation 
A c /Ap — 0.003(n — 1.75) 3 . Crosses: Cox et al's experiment [5], 
the mean critical value at which the central bubble is ejected 
when its area decreases (A**). The dashed line is the quadratic 
fit from [5], as described in text. Dotted, the theoretical curve 
where the dry flower cluster would undergo symmetry reduc- 
tion [23]. 

Figure [16] compares the values of Au m and the exper- 
imental A** [5] as functions of n. Our values are much 
closer to the experimental data than those of the symme- 
try reduction criterion. The calculated values of A^ fit 
very well to a cubic power law: 

A nm = 0.003(n - 1.75) 3 ^ p . (43) 

This result disagrees with the quadratic law in |5J. This 
quadratic law was found by adjusting the central area so 
that the energy of the symmetric cluster and that of the 
ejected configuration be equal, as calculated by Surface 
Evolver [30] . 

In figure [TBI the small, but systematic, difference be- 
tween the experimental ejection area and An m may be 
explained by considering that, experimentally, the size of 
the bubbles can not be exactly the same, thereby reduc- 
ing the shortest film length and triggering Tl a bit earlier 
than in an exactly symmetric model. 

Figure [17] shows the predicted liquid fraction at Au m as 
a function of n. Here, the liquid fraction is defined as the 
ratio of the liquid to total (enclosed gas+liquid) area|f|. At 
constant pressure, the liquid volume fraction has to vary. 
The experiments appear to have been performed at about 
4% liquid fraction^. 



4 This $ may differ from global liquid fractions in not in- 
cluding any gas or area outside the cluster. 

5 This is an area fraction in the 2D model. Specifying the 
physical liquid volume or fraction is more delicate, if relevant 
at all, due to the finite thickness in the third dimension and 
the liquid borders along the plate(s). 
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Fig. 17. Values of the liquid fraction obtained at the limit 
decoration versus the petal number n f figure [16)) . 



7 Conclusions 

To summarise, in a 2D equilibrated foam, star-triangle 
equivalence permits to replace dry three-fold vertices by 
three-sided cells and vice versa. The general structure of 
a decorated star was obtained by first defining a two real 
parameters reference 3-cell, then applying a one complex 
parameter Moebius transformation. This way allowed us 
to give to the Plateau border and 3-sided bubble problems 
a unified description. We calculated all the geometrical 
and physical quantities of interest: area, perimeter and 
energy. Beside Plateau borders and standard three-sided 
bubbles, the method also applies to triple partitions of the 
disk and the three-bubble cluster. 

The difference between the decorated (triangle) and 
bare (star) vertex is quantified by the relative excess en- 
ergy: _ 

the energy difference normalised by the square root of the 
triangle area. This dimensionlcss quantity, depending on 
only three parameters, serves as a shape discriminant. At 
given area, the most symmetric configuration minimises 

Fixing the equilibrated star, the relative excess energy 
of 3-cells is affine in A with a good approximation: 

e (A)~e° + bA, 

where b depends on the form of the decorated films. The 
slope b is, up to a positive factor, the mean square curva- 
ture of the star (K^ ta ); being zero in the 3-fold symmetric 
case, b measures the deviation from perfect symmetry. 

As an application to the flower problem, we calculated, 
for various petal numbers and including liquid at the ver- 
tices, the limit value of the central area when two Plateau 
borders get into contact. Tuning the liquid pressure to fit 
the experimental data, we obtained a nice agreement of 
the calculated threshold with the experiment. We predict 
an experimental liquid fraction of about 4%. 

The calculated threshold is weakly sensitive to the liq- 
uid pressure. Very likely, the small systematic difference 
between the calculated and experimental data comes from 



a slight discrepancy of the experimental bubble size with 
exact n-fold symmetry. The calculated limit central area, 
as a function of the number of petals n, fits very well a 
cubic law. 
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Appendixes 

A Moebius transformations 

A Moebius transformation or homography, z = f(z), is 
an automorphism of the augmented complex plane (C* = 
C U {oo}), / : C* -> C* [31132] . defined by: 



/(*) 



a z + b 

cz + d 



with ad — be = 1 . 



(44) 



A Moebius transformation is a combination of transla- 
tions, rotations, dilatation and inversions 3 lj . 

Beside preserving angles (conformality), these trans- 
formations map circles to circles; in particular any circle 
passing through the inversion point, v* = —d/c, is sent to 
a straight line (circle through infinity). 

The circle of parameters (zq = 0, r) £ C x R (centre 
and radius) is mapped by (144|) to the circle of parameters: 



(zo,r) 



acr 



bd \ad — bc\ r 



| c |2 r 2 _ | d |2> | r 2 | c |2 _ | d |S 



(45) 



where x stands for the complex conjugate of x. 

Similarly, / maps the straight line of parametric equa- 
tion z = zq + rat, t G U C K, with zo = , m G C, to the 
circle of parameters: 



2alm(racd) + icm \ad — b( 



2 c Im(c ra d) 



1 2 Im(cm d)\ 



(46) 



When z ^ 0, equations ([45]) and (f46|) keep their validity 
with the substitution: b — > b + qzq and d — > d + czq. 



B 3-cell perimeters 

Let us consider an edge defined by its parametrisation 
f(u), u G [uq, ui]. In the transformed plane, the length of 
tp = f s o (p, image of <p by / s , is 



I = I |^(m)| dtt 



(i-N 2 ) 



u \l-S(p(u)\ 



ip(u)\ du = 
y|#tO|du. (47) 
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B.l Internal perimeter 

For an internal film <^ ta (ii), substituting equation ((8j) in 
([4"T[) gives, after simplification: 



If* = (l-\s\ 2 ) 



-. - 



dit 



rf a (s)Qf*(s,q), (48) 



l Q l-2u|s|cos(6» s -^-)+w 2 |s| 2 

t; 

where u)f a (s,q) is given in ([22 



B.2 Triangle perimeter 

For a triangle edge ipffe), (|4"7|l can be reduced to an 
integral of the form: 



rt ri _ g(l-M 2 ) ( l 
j 2 I sin/3 I 



(3du 



i aj + Oj cos uf3 + a J 2 sin u(3 



(49) 



where the constants in the integral are: 
\s\q 



V3 



cos( 



% - e s ){l - VScot/3) + 



1 1 2 (5 + 2 cos(2/3 + 7r/3)), (50) 



12 sin 2 /3 



-.3 — 



\*Q 
sin/? 

\s\q 



») + » m ■ 7, ' ( 51 ) 



2a/3 sin /3 



sin 







sin/3 



they verify 



K) 2 - (a]) 2 - (a J 2 ) 2 - q2{ \J S][2 f > . 
V 01 V 11 V 2/ 4 (f ") 2 sin 2 /3 



(52) 



(53) 



Let A be the square root of the previous expression. Then 
the solution of the integral (|4"9"]h verifying (fB"3"]) . is 33J : 



rtri = g(l-|s| 2 ) 
3 A\smp\ 



arctan 



(a 3 Q — a{) tan(|) + a J 2 
A 



+0 



Taking, for A, the root of the second member of (|53|) gives 
lf = ff (s,<7,/?)^f (s,q,l3), (54) 



if = 2 arctan(2ff sin (3 



(55) 



Equations ([541 [55)1 also define the angle o)j n subtended by 
the edge (pj 1 . 

C Series coefficients of ep 

The coefficients in the expressions ([3"5"|) . ([3^)1 and (|3T[) are 



if a = %/3, U 



3-y/3 



tri 



/8 



sin/3 ' 



7- tri 

^3 — 



^2 = 



+ 5/3 - 2%/3sin(2/3 + f ) - 2/3sin(2/3 - |; 



1 



3/3 



4sin 3 (/3) 
2 



-4 4 



4sin/3 \sin/3 cos(/3 + 7r/6) / 
1 f 4( V3 -0) -(9 + 4V3 /3) cot /3 



esc 2 /3 + 

+ - /3csc 4 /3 -=(1 + a/3 cot /3) . 

8 M v 4V3 

The zeroth and second order (in q) coefficients of L tri and 
i sta are zero at q = 0. In the same way Aq, Ai and A3 
are zero. 
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